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1 Introduction 

Two of the most fascinating problems in astrophysics and cosmology today 
are the origin and nature of the dark matter (DM) and the sources and 
mechanisms behind the so-called gamma-ray bursts (GRB). The presence of 
non- luminous, gravitationally interacting matter was inferred in the earlier 
part of this century when the dynamics of groups of galaxies did not agree 
with the predictions based on the observed amount of luminous matter. Ever 
since, physicists have tried to build models that can explain both the nature 
and the origin of the dark matter, in terms of either hot dark matter (HDM) 
or cold dark matter (CDM), referring to their internal velocity distribution. 
Recent analysis [1] seems to favour a mixture of these two categories, at least 
when compared to the observed cosmic microwave background radiation 
(CMBR) and the standard theory for the large scale structure formation 
in the Universe (arising from gravitational enhancement of small density 
perturbations in the early Universe). 

GRBs were first observed by the VELA satellites in the late 1960s [2]. 
These bursts, isotropically distributed across the sky, are assumed to orig- 
inate from cosmological distances, and they radiate up to 10 54 ergs in a 
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time duration of 0.01 — 1000 s. They are the most luminous sources in the 
Universe. Several models have been proposed for the origin of the bursts, 
and as of today, none has been singled out as most favoured. So-called fire- 
ball models (see [3] for an overview) where a relativistic expanding shell, 
originating from the source, interacts with the intergalactic matter is very 
popular since they can explain the characteristics of the radiation. However, 
the actual source, the "inner engine" , of the energy release is not well un- 
derstood. The most common explanation includes a merger of two compact 
objects of some kind (neutron stars, black holes etc.). 

Our work is built on one, basic principle: to use known physics to explain 
new phenomena. Our starting point has been the widely accepted theory 
that in the early Universe, t < 10 -5 s after the Big Bang, all matter was in 
the form of quark-gluon plasma (QGP). At t « 10~ 5 s, a phase transition 
occurred, which confined the free quarks into hadrons. This quark-hadron 
transition is assumed to be of first order, something that several lattice 
calculations seem to indicate [4]. In 1984 Witten [5] suggested the possibility 
that regions of the high temperature phase QGP would remain and be stable 
after the transition. If this is correct, then a very exciting possibility to 
explain both the origin of the baryonic dark matter and the sources for 
GRBs opens up. In our model, we identify the sources of GRBs with the 
baryonic dark matter. 

The fundamental assumption that underlies our model is that the QGP 
is the absolute ground state of QCD. This idea is not new, in fact it emerged 
in the early 1970s. An admixture of u, d and s quarks is likely to have a 
lower energy per baryon number than the proton. This makes it possible 
for the QGP to be stable, even on cosmological scales. 

To summarize, we make two major assumptions in our model: 

1. QGP is the ground state of QCD, at least for massive objects. 

2. The quark-hadron transition was of first order. 

2 The model 

The starting point for our model is the assumption that only a very small 
part of the primordial QGP actually hadronised into ordinary matter. The 
phase transition allowed QGP objects survive the transition, or "quark 
nuggets" as other authors call them, fractally distributed in size and oc- 
curring as inhomogeneities in the surrounding mixture of hadronic matter 
and vacuum. Smaller QGP objects, with baryon number A < 10 44 disap- 
peared due to neutron and proton condensation [6], made possible by the 
high temperature, while larger objects survived. 

Since the average density of one of these quark objects should be around 
10 13 g/cm 3 , and the radius of an object with 10 44 quarks is on the or- 
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der of a centimetre, gravity becomes important when deciding under what 
circumstances the QGP objects are stable with respect to gravity and to 
the internal degeneracy pressure occurring due to the Pauli principle. The 
mass-radius relationship of a large spherical QGP bag (A > 10 50 ) can be 
calculated with the Tolman-Oppenheimer-Volkoff equations [7], derived in 
general relativity: 



dp [e(r) + p(r)\[m{r) + 4irr 3 p(r)] 
dr r[r — 2m(r)] 

dm{r) 



(1) 



= 47rr 2 e(r) m(r) = Air f e(r')r' 2 dr' (2) 
dr J 

p(r = 0) = p c (3) 



p(r = R) = 0. (4) 

Solving these equations (with c = G = 1) requires an equation of state 
for the QGP. Our result is based on treating the QGP as a relativistic fermi 
gas with zero chemical potential and no interactions. The phase boundary is 
expressed in the spirit of the MIT bag model with an external bag pressure 
characterised by the bag constant B: 

e(r) = 3p(r) + AB. (5) 

The kinetic part of the pressure is (with three quark flavours /): 

^^ + £(1^ + 1x^ + ^4 (6) 

When integrating the differential equation (1), a mass-radius relationship 
emerges, as shown in Figure 1. Using these calculations for deciding the 
upper mass limit, QGP objects surviving the quark-hadron transition and 
hadron evaporation must lie in the range 1CT 17 M Q < iiiqgp < 2 M with 
radius 10~ 5 km <tqgp < 10 km, or 10 44 < A< 10 58 . M is the solar mass 
2-10 33 g. 

The size of the surviving QGP objects is characterised by the nucleation 
distance in the phase transition. Ref. [9] suggests a nucleation distance on 
the order of a cm, while others [10] allow distances on the order of a km. The 
assumption that the transition occurred through detonation mechanisms 
[11] has not been disproved, making it impossible to exclude km-sized QGP 
remnants. The horizon size at the transition is on the order of km with an 
enclosed mass of roughly M . It is an interesting fact that the mass energy 
enclosed in a solar-mass QGP object is roughly the same as the gamma-ray 
energy in a GRB at cosmological distances, z = O(l). 
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Figure 1: The stability relations (full lines) between the mass and the radius 
of an spherical QGP object for different values of the external pressure B. 
B 1 / 4 = 180 MeV (curve 1.), 150 MeV (2.), 120 MeV (3.), 100 MeV (4.), 
75 MeV (5.) and 50 MeV (6.). For clarity, only the most relevant segments 
of the full lines are shown. The hatched line shows the criterion for collapse 
into a black hole, as given by the Schwarzschild radius. Results for other 
parameter values can be found in [8j. 

If the QGP objects survived the phase transition and the evaporation 
they should be stable for a cosmological period of time, making up for at 
least the baryonic part of the cold dark matter. The only way they could 
hadronise or decay is if an external energy is added. One such circumstance 
is if two of these objects merge and add gravitational energy to the system. 
This is analogous to how a supernova can make elements heavier than iron. 
When a QGP hadronise different radiative mechanisms [12] might contribute 
and gamma-ray emission is one of them. This is why we suggest that gamma- 
ray bursts come about when parts of the QGP objects hadronise. One should 
notice though that the mean free path for a gamma-ray in a quark-matter 
environment is roughly 100 fm, which poses some difficulties in explaining 
how such a vast amount of energy, 10 54 ergs, is radiated during a merger. 
This is under further investigation, but one should notice that this problem 
is shared by all proposed mechanisms for mergers as sources for GRBs. Our 
scenario includes all the possible mechanisms for gamma-ray production 
stated in the conventional scenarios and, in addition, includes the possibility 
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of direct gamma-ray radiation through hadronisation. 

One crucial point is of course the number density of QGP objects that 
merge. We know from observations that the average galaxy hosts approx- 
imately 1CT 4 — 1CT 6 bursts/yr and that must be made to fit our initial 
distribution of solar mass sized QGP objects. An important fact is also that 
all bursts seem to originate at cosmological distances. These two observa- 
tional facts suggest that mergers were more frequent in the distant past and 
that the fraction of QGP objects (or at least of binary systems) that lies in 
the M range is fairly small. 

We can conclude that there are no observational data that contradict 
the idea of QGP objects as the sources for GRBs. Another consequence is 
that of "beaming" of the gamma rays, since the hadronisation in a merger 
probably occurs through a bridge between the two objects. This also lowers 
the required energy alleviating the problem mentioned above. This bridge 
forms when the perturbing gravitational potential from the partner in the 
merger is large enough to make it possible for the QGP to hadronise. This 
means that the largest perturbation occurs along the symmetry line between 
the centres of the two objects. 

3 Particle physics constraints 

One of our main assumptions is that the QGP is the true ground state of 
QCD. This means that the energy content per quark is lower in a QGP 
than in a nucleon. Such ideas began to flourish in the early 1970s with a 
pioneering work by Bodmer [13]. Since then several calculations supporting 
this assumption have been made. Farhi and Jaffe [14] explored the properties 
of the so-called strange matter and showed that within certain parameter 
ranges, the energy per baryon in strange matter is lower than in the proton. 
Witten [5] showed that the presence of a s quark in a u, d quark mixture 
lowers the energy per quark by a factor of 0.89 as compared to non-strange 
matter. Other analyses of the stability of strange matter were performed 
by Bjorken and McLerran [15], Chin and Kerman [17] and De Rujula and 
Glashow [18]. They suggested that since the mass of the s quark is lower 
than the Fermi energy of a u or d quark in a multiquark system, the opening 
of a new flavour degree of freedom could lower the Fermi energy of the 
system. This behaviour can be seen in the simulation shown in Figure 2. 
Other simulations have shown similar results but no such stable multiquark 
objects have been detected. The existence of such a quark matter phase is 
hence just a theoretical prediction. The reason why atomic nuclei do not 
decay into a QGP is that the conversion of a u or a d quark into an s quark 
is a weak process with a very low probability. The transition of 56 Fe into 
strange matter containing 168 u, d and s quarks in equal proportions is a 
56th order weak process. Hence, the probability of decay of ordinary matter 
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Figure 2: This plot shows the dependence of the energy "per baryon" of a 
giant bag on the baryon number for = 145 MeV [19]. The matter in the 
bag is SU(3) flavour symmetric and electrically neutral. 

into strange matter is exceedingly small. It should be noted that it is not 
obvious that models of confinement, such as the MIT bag model [16], can be 
used for bags of arbitrary sizes or that the bag constant must be independent 
of the number of confined quarks. 

The details of the quark-hadron transition are not very well known. 
Large experiments (either running or planned) will try to detect the phase 
transition. The consensus among both cosmologists and particle physicists 
is that it occurs at a temperature of ~ 150 MeV and is of first order. Lattice 
calculations seem to indicate a first order transition at that temperature but 
it is fair to say that the results are inconclusive, due to tremendous difficul- 
ties of simulating a fermion system with a non-zero chemical potential. This 
applies both to the type of the transition, the temperature and the mean 
nucleation distance. 

The possible survival of QGP objects has been discussed for some time. 
Some authors have argued [20] that large QGP objects would boil when 
hadronic gas is formed in their interior due to the superheating they are ex- 
periencing. However, the importance of this mechanism has been questioned 
[21], since the nucleation of hadronic bubbles in the interior is a rather slow 
process compared to surface evaporation and cooling. 
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4 Astrophysical constraints 



Two different observations will be discussed here. One of these is the CMBR 
which comes from the last scattering surface, ~ 10 5 years after the Big Bang. 
During the last nine years, the COBE satellite has measured the CMBR all 
over the sky. The data describe an almost isotropic, black-body radiation 
at a temperature of T = 2.725 ± 0.005 K [22]. The small anisotropies in the 
radiation distribution, ~ 30 ± 5 //K, is assumed to originate from small 
density perturbations present at the time of the last scattering, since the 
anisotropy is related to fluctuations in the density distribution, 

T p v ' 

If our scenario of surviving large QGP objects is correct, then one might 
think that this is not compatible with the COBE data due to the huge 
"lumpiness". This is not the case since only non-neutral matter interacts 
with the photons present at the time of last scattering. Since the QGP 
objects are charge-neutral made up of u, d and s quarks and a small amount 
of electrons, the possible imprint left on the CMBR by the QGP objects is 
most likely unobservable. We would like to add, however, that this is not a 
settled issue and a thorough investigation of the exact interaction between 
the CMBR and the QGP objects remains to be made. 

The other major observation we want to address is the abundances of 
light elements. The standard theory ("SBBN") describing the formation of 
these elements states that they appeared when the Universe was between 
10~ 4 s and a few minutes old. Since the QGP objects in our model are in 
their ground state, these did not partake in the formation of light elements. 
The presence of large inhomogeneities in the nucleosynthesis era affects the 
nucleosynthesis. Various authors [23-26] have addressed inhomogeneous Big 
Bang nucleosynthesis and recent results [27] indicate that the presence of 
inhomogeneities on the metre to kilometre scale actually reduces the discrep- 
ancy between the observed abundances of 4 He and D and the prediction of 
SBBN. Ref. [27] states that the origin of these inhomogeneities could very 
well be the quark-hadron transition. 



5 Conclusions and outlook 

One can conclude that no observations rule out the possible presence of 
relic QGP objects, acting as both the baryonic cold dark matter and as 
the source for gamma-ray bursts. However, one should also say that no 
direct observations support this idea. The two assumptions we have made 
regarding the stability of QGP and the nature of the phase transition have 
no experimental evidence either. There is no doubt that our model is to 
be considered as both rudimentary and speculative. Much work remains to 
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be done, where perhaps the most important one is pure particle physics; to 
establish whether or not the QGP is the ground state of QCD. The results of 
future experiments at RHIC and LHC will hopefully confirm the existence 
of a deconfined phase, the QGP. On the theoretical side, a detailed analysis 
of the split-up of the QGP phase in the quark-hadron transition is desireble. 

On the observational side there is certainly an enormous amount of in- 
formation in the characteristics of the CMBR remaining to be discovered, 
which improved data taking and analysis should be able to extract. For 
this, a more detailed analysis of the possible interaction between the QGP 
objects and the CMBR is necessary. 

This work has been done in collaboration with Sverker Fredriksson and 
Johan Hansson at the Department of Physics, Lulea University of Technol- 
ogy, Svante Ekelin at the Department of Mathematics, Royal Institute of 
Technology, Stockholm and Argyris Nicolaidis at the Department of Theo- 
retical Physics, Aristotle University of Thessaloniki. 
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